Eigenvalue Distribution In The Self-Dual Non-Hermitian Ensemble 
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We consider an ensemble of self-dual matrices with arbitrary complex entries. This ensemble is 
closely related to a previously defined ensemble of anti-symmetric matrices with arbitrary complex 
entries. We study the two-level correlation functions numerically. Although no evidence of non- 
monotonicity is found in the real space correlation function, a definite shoulder is found. On the 
analytical side, we discuss the relationship between this ensemble and the = 4 two-dimensional 
one-component plasma, and also argue that this ensemble, combined with other ensembles, exhausts 
the possible universality classes in non-hermitian random matrix theory. This argument is based 
on combining the method of hermitization of Feinberg and Zee with Zirnbauer's classification of 
ensembles in terms of symmetric spaces. 



PACS numbers: 05.45+b, 73.20.Dx, 03.65. Ca 
I. INTRODUCTION AND CLASSIFICATION 

There are ten known universality classes of hermitian 
random matrices. Dyson jlj proposed the existence of 
three symmetry classes, depending on spin and the exis- 
tence of time reversal symmetry. These give the three 
classes known as Gaussian Unitary, Orthogonal, and 
Symplectic (GUE, GOE, GSE). Another three ensem- 
bles are the chiral Gaussian ensembles (chGUE, chGOE, 
chGSE) . These ensembles are of relevance to low en- 
ergy QCD. Altland and Zirnbauer introduced four more 
ensembles which can appear in superconducting systems 
||. Finally, Zirnbauer demonstrated a relationship be- 
tween the different classes of random matrix theory and 
symmetric spaces, and from this argued that the ten dis- 
tinct known universality classes exhausted all possible 
universality classes Q . 

In this section we discuss various universality classes of 
non-Hcrmitian random matrices, including the ensemble 
of arbitrary self-dual matrices, the subject of this paper. 
We mention the concept of weak non-Hermiticity, but do 
not consider it further in this paper. We argue that the 
various classes of non-Hermitian matrices, the self-dual 
ensemble and four others, exhaust all possible universal- 
ity classes. Finally, possible applications of the self-dual 
ensemble are dealt with, including relations with the one- 
component plasma. In section II, we further discuss the 
relationship with the one-component plasma. In section 
III, numerical results for the self-dual ensemble are dis- 
cussed, in particular the eigenvalue density as a function 
of radius and the two-eigenvalue correlation functions. 

Several ensembles of non-Hermitian random matrices 
are common in the literature. Ginibre Q introduced 
three classes of such matrices. One is an ensemble of ma- 
trices with arbitrary complex elements, one an ensemble 
with arbitrary real elements, and the third an ensemble 
with arbitrary real quaternion elements. Another ensem- 
ble of non-Hermitian matrices is an ensemble of complex, 



symmetric matrices || . This ensemble arises particularly 
in problems of open quantum systems. This gives a total 
of four known universality classes. 

For each of these ensembles, there exists a weakly non- 
Hermitian version of that ensemble. This idea of weak 
non-Hermiticity was introduced by Fyodorov et al. ||. In 
this case the anti-Hermitian part of the matrix is small; 
we only consider strongly non-Hermitian matrices in the 
present paper and do not consider weakly non-Hermitian 
matrices, even though they are the most relevant for scat- 
tering problems. 

The strongly non-Hermitian ensembles can be obtained 
from a general three parameter family of non-Hermitian 
matrices introduced by Fyodorov et al. [Q. This family 
includes parameters measuring the strength of the real 
and imaginary, symmetric and anti-symmetric parts of 
the matrix. By adjusting the parameters, one can obtain 
various ensembles. One possiblity, which does not appear 
to have been considered much, is an ensemble of anti- 
symmetric matrices with arbitrary complex elements. 

Now, let us show that this ensemble is equivalent to 
an ensemble of self-dual matrices with arbitrary complex 
elements; this is the ensemble considered in this paper. 
Let A be an arbitrary anti-symmetric matrix. Let Z be 
the matrix given by 
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(i) 



Then, Z T = -Z and Z 2 = -1. Let M = ZA. It is 
trival to verify that ZM T Z = —M. So, M is self-dual. 
The advantage of using self-dual matrices instead of anti- 
symmetric matrices is that self-dual matrices have pairs 
of equal eigenvalues while anti-symmetric matrices have 
pairs of opposite eigenvalues; this makes the correlation 
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functions clearer. When choosing matrices from the en- 
semble, we will use Gaussian weight 

e _iTr(M+M) ( 2 ) 

Given these five classes, the 3 ensembles of Ginibre as 
well as the ensembles of symmetric non-Hcrmitian and 
self-dual non-Hermitian, let us ask whether all possible 
universality classes of strongly non-Hermitian random 
matrices have been found. Feinberg and Zee introduced 
the method of hermitization for non-Hermitian matrices 
H . A similar technique was used by Efetov jl0| . The ba- 
sic idea is to take a non-Hermitian matrix M — E, where 
E is a complex number, and form the Hermitian matrix 

jr = ( M h - E R M a - iE X \ 

\Ml+iE! -M h + E R ) {6) 

where Mh is the Hermitian component of M and M a is 
the anti-Hermitian component of M and En and Ei are 
the real and imaginary components of E. Equivalently, 
one can form the Hermitian matrix 

h= {m^-e o ) W 

From the zero eigenvalues of H, one may extract the 
zero eigenvalues of M — E. So, to each universality class 
of non-Hermitian random matrices, there corresponds a 
universality class of Hermitian random matrices. 

If we hermitize the three non-Hermitian ensembles in- 
troduced by Ginibre, we obtain the three chiral ensem- 
bles (chGUE, chGOE, chGSE). The relation with the chi- 
ral ensembles is most clear using equation (||), instead 
of equation (Q). If we hermitize the ensemble of sym- 
metric, complex matrices we obtain the ensemble with 
symmetry class CI, according to the nomenclature of Al- 
tland and Zirnbauer. If we hermitize the ensemble of 
self-dual complex matrices, we obtain the ensemble with 
symmetry class DHL Here the relation with the Hermi- 
tian ensembles is most clear using equation (^). The 
other five classes of hermitian random matrices cannot 
be obtained by hermitizing a non-Hermitian ensemble: 
the GOE, GUE, and GSE classes lack the needed block 
structure, while the C and D ensembles lack the symme- 
try that relates the elements in the upper left and lower 
right blocks. This suggests that all possible universality 
classes of non-Hermitian matrices have been obtained. 

One possible interest in the ensemble of self-dual com- 
plex matrices would be experimental, such as in open- 
systems with spin orbit scattering. Another interest is 
theoretical, considering the relationship of this ensemble 
to the (3 = 4 one-component plasma in two-dimensions. 
Although the level distribution in the ensemble differs 
from the distribution of charges in the plasma, there are 
some close relations between the two, discussed more in 
the next section. 



It is known that the ensemble of matrices with arbi- 
trary complex elements is equivalent to the (3 = 2 plasma. 
The correlation function of the (3 = 2 system is mono- 
tonic, with Gaussian decay. From perturbative calcula- 
tions pH , it has been suggested that, for (3 > 2, the two- 
level correlation function becomes non-monotonic, indi- 
cating the appearance of short-range order. This makes 
it very interesting to examine the correlation function of 
the ensemble of self-dual matrices, although no signifi- 
cant sign of any non-monotonicity is found here in the 
numerical calculations. 

Numerical calculations on the one-component plasma 
|l~2| suggest that there is a phase transition at j3 » 144; 
so, any order that exists for (3 = A must be short 
range. An exact study for finite number of particles |l~3j | 
showed non-monotonicity of the correlation functions for 
(3 = 4, 6. Even for (3 = 4 there is a definite peak in the 
correlation function. 

II. (3 = 4 ONE-COMPONENT PLASMA 

Consider a system of N particles, located at positions 
Zi, with partition function 

r N 

dz^We-^We^ ^-^ (5) 

i=l i<j 

This defines the two-dimensional one-component plasma. 
For (3 = 4, there exists some relation between this system 
and the ensemble considered here. 

First, the density of the plasma, p, is equal to 
where the density is measured in charges per unit area. 
The plasma has constant charge density p in a disc about 
the origin, and vanishing charge density outside. The 
self-dual ensemble has the same charge density, as found 
numerically in the next section, and as can be shown with 
a replica or SUSY technique. 

Second, there exists a relationship between the joint 
probability distribution of the eigenvalues of M and the 
probability distribution of charges in the one-component 
plasma. The j.p.d. of the eigenvalues of M is different 
from the charge distribution in the plasma, but we will 
argue that for widely separated eigenvalues the j.p.d. of 
the eigenvalues behaves the same as the probability dis- 
tribution of the charges. 

Let M be a matrix within the ensemble of self-dual, 
complex matrices. We can write M as M = AAA -1 , 
where A is a diagonal matrix of eigenvalues of M. The 
eigenvalues of A exist in pairs, with [A, Z] = 0. The 
requirement that M be self-dual is equivalent to the re- 
quirement that X T Z = ZX- 1 and Z(X- X ) T = XZ; 
if this constraint on X holds it is easy to verify that 
ZM T Z = -M. 

If we were to impose the additional constraint on X 
that X be unitary, then we would find that X must be 
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an element on the symplectic group. In this case, with X 
in the symplectic group, the matrix M must be normal, 
such that [M, M<] =0. In this case the distribution of 
eigenvalues of M exactly matches the charge distribution 
in the (3 = 4 plasma. 

In the general case, M is not normal and X is not 
unitary, and the distribution of eigenvalues of M will be 
different from the charge distribution of the plasma. Still, 
consider a situation in which we fix A and integrate over 
X, with Gaussian weight e -2 Tr ( M M ). This is how one 
obtains the j.p.d. of the eigenvalues. 

The measure [dM] on matrices M is equivalent to the 
measure [dAj][dX] Y[ — -\?| 8 . The j-P-d. of the eigen- 

i<j 

values is defined by 



in* 



A, I 8 



[dX}e 



-iTrtM+M) 



(6) 



-iTr(MtM) 



with M = XAX^ 1 . The Gaussian weight, e" 
will depend on X. It will be greatest when X is chosen 
to be symplectic, so that M is normal. If the eigenval- 
ues of A are well separated, then the exponential in the 
Gaussian weight will be large, and we can evaluate the 
integral by a saddle point method: we will restrict our 
attention to a saddle point manifold of matrices M which 
are normal, as well as weak fluctuations away from this 
saddle point manifold. If we parametrize the fluctuations 
away from the saddle point manifold and then treat these 
fluctuations in a Gaussian approximation, valid when the 
eigenvalue separation is large, we obtain that the j.p.d. 
for the self-dual ensemble is equal to, in this particular 
approximation, 



TV 



II e " |Zl1 11(1* -Zj^Hdzidz 



(7) 



l<3 



up to a constant factors. This is, of course, the same 
as the probability distribution of the charges in the one- 
component plasma at (3 = 4. 

In general, we expect that for well separated eigen- 
values, the level repulsion in the self-dual ensemble will 
match the charge repulsion in the plasma; it is only the 
short distance interaction that will be different. Further, 
it may be shown explicitly by calculations on small ma- 
trices that the short distance interactions in the j.p.d. for 
the self-dual ensemble cannot be written as a product of 
two-body terms. 

Given these similarities, one might hope that the corre- 
lation functions of the self-dual ensemble will shed some 
light on correlations within the plasma. In the next sec- 
tion, we discuss a numerical investigation of the self-dual 
ensemble. 



III. NUMERICS 

Mathematica was used to generate 4940 600-by-600 
self-dual matrices. The matrices were chosen with Gaus- 



sian weight e 



-iTr(M+M) 



as in equation 



The matrices 



have 300 pairs of eigenvalues. A picture of these eigen- 
values for a typical matrix is shown in figure 1. 
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FIG. 1. Plot of eigenvalues for a typical 600-by-600 matrix 

The eigenvalue density as a function of radius is shown 
in figure 2. The density obeys the circular law dJT^]: it 
is nonvanishing and roughly constant within a disc, and 
vanishing outside. For the /3 = 4 one component plasma, 
with a confining potential e~ zz (see equation (||), the ex- 
pected density of particles per unit area, from the circu- 
lar law, is The single particle eigenvalue density ob- 
served numerically for the self-dual matrices agrees with 
this result; note that since eigenvalues come in pairs, then 
we expect e~ zz to be the confining potential that corre- 
sponds to the weight of equation (||) as each eigenvalue 
in the pair contributes a factor of e~ . 
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FIG. 2. Average eigenvalue density as a function of radius 

One interesting feature of figure 2 is that the eigenvalue 
density near the edge rises before dropping. It is not clear 
why this happens. 

The two-level correlation function is shown in figures 
3 and 4. In figure 3, we look at all eigenvalues within a 
distance of 6 or less from the origin, and plot the prob- 
ability to find another eigenvalue at given distance from 
the first eigenvalue. In figure 4, to reduce effects due to 
the finite size of the matrix M, we require that the first 
eigenvalue lie within a distance of 3.5 or less from the ori- 
gin. No significant differences are found between figure 3 
and figure 4, indicating that the effects due to the finite 
size of M are small even in figure 3. 



eigenvalue density near the edge, as shown in figure 2, 
and has no deep meaning. 

Looking at figures 3 and 4, there is a definite "shoul- 
der" at a distance of slightly less than 3. There is no def- 
inite sign of any non-monotonicity; certainly, if there is 
any peak in the correlation function near the shoulder, it 
is much smaller than the peak found in the j3 = 4 plasma 
p3[ . As a quick estimate of the expected spacing between 
levels, assume that the levels formed a perfect hexago- 
nal lattice, so that they are very ordered, and packed as 
closely as possible. In this case, if the levels have a den- 
sity of , then the closest spacing between levels is , 
which is approximately 2.7. For other arrangements of 
levels, the spacing will be slightly less. This length agrees 
quite well with the size of the shoulder. So, the shoulder- 
length matches reasonably with the length scale expected 
from the particle spacing. 



IV. CONCLUSION 

In conclusion, we have considered an ensemble of 
strongly non-Hermitian, self-dual matrices. The two- 
level correlation function of this ensemble is particularly 
interesting, although the hoped for non-monotonicity has 
not emerged. It seems that all possible universality 
classes of non-Hcrmitian matrices are now known. 
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FIG. 3. Average two level correlation function. See text. 
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FIG. 4. Average two level correlation function. See text. 

One can see finite size effects in both figure 3 and 4 
for large distances. The correlation function rises for dis- 
tances of around 20. This is simply due to the rise in 
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